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Abstract 

Estuarine Batagur are among the most critically endangered chelonian species. We assess the taxonomic status of the 
recently discovered Cambodian relic population of Batagur by phylogenetic analyses of three mitochondrial (2096 bp) 
and three nuclear DNA fragments (1909 bp) using sequences from all other Batagur species and selected allied 
geoemydids. Furthermore, we calculated haplotype networks of the mitochondrial cytochrome b gene for Cambodian 
terrapins, B. affinis, B. baska, and B. kachuga and compare external morphology of estuarine Batagur populations. 
Genetically, Cambodian Batagur are closely related with, but distinct from B. ajfinis from Sumatra and the west coast of 
the Malay Peninsula. Morphologically, Cambodian Batagur resemble the distinctive B. affinis populations from the 
eastern Malay Peninsula that were not available for genetic study. We suggest that the Batagur populations from the 
eastern Malay Peninsula and Cambodia represent a new subspecies of B. affinis that once was distributed in estuaries 
surrounding the Gulf of Thailand {Batagur affinis edwardmolli subsp. nov.). Its patchy extant distribution is most 
probably the result of large-scale habitat alteration and century-long overexploitation. In addition, our phylogenetic 
analyses suggest repeated switches between riverine and estuarine habitats during the evolution of the extant Batagur 
species. 

Key words: Southeast Asia, South Asia, Batagur affinis affinis, Batagur affinis edwardmolli subsp. nov., Batagur baska, 
Batagur kachuga, endangered species 


Introduction 

Batagur baska (Gray, 1830) is one of the world’s most critically endangered terrapins. Its range was recently 
restricted to a region extending from coastal north-easternmost India and adjacent Bangladesh southwards to 
at least the Ayeyarwady and Bago estuaries in Myanmar (Praschag et al. 2007, 2008). Populations from the 
Malay Peninsula and Sumatra, traditionally treated as conspecific, turned out to represent the distinct species 
B. affinis (Cantor, 1847). Both B. affinis and B. baska are large terrapin species, reaching a maximum shell 
length of approximately 60 cm (Moll 1980; Ernst et al. 2000). Their distribution is more or less confined to 
brackish water; they occur in estuaries, mangrove belts and inshore beds of marine vegetation (Kalyar et al. 
2007). Phylogenetically, B. affinis and B. baska together are sister to a riverine inland species, B. kachuga 
(Gray, 1831) from northern India (Praschag et al. 2007), that was placed for a long time in the genus Kachuga 
(Ee et al. 2007; Praschag et al. 2007). 

Batagur affinis and B. baska declined dramatically throughout their ranges as a result of overharvesting of 
adults and eggs coupled with habitat degradation (Das 1997; Moll 1997; Kalyar et al. 2007; Platt et al. 2008). 
Estuarine Batagur once also occurred in southern Vietnam and Cambodia, but until now, nothing has been 
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known about their taxonomic identity. According to Bourret (1941), ‘5. baska" was common in the late 19* 
Century in Cochinchina (southern Vietnam) and Cambodia, and eggs were collected along the Prek Tap 
Cheang River, now known as Sre Ambel, a coastal river in the province of Koh Kong, Cambodia (Platt et al. 
2003). The latter authors added that ‘R. baska' also occurred in the early 20* Century in the Tonic Sap, a large 
freshwater lake in the inland of Cambodia. While Platt et al. (2003) failed to confirm the species for the Tonle 
Sap, an unexpected habitat for an estuarine species, they were able to provide evidence for the existence of a 
relic population in the Sre Ambel River system. In the present study we compare specimens of this Sre Ambel 
population genetically with B. baska and B. affmis using sequence data of three mitochondrial (2096 bp) and 
three nuclear DNA fragments (1909 bp). Further, we compare external morphology and natural history of 
estuarine Batagur populations, correlate these data with observed genetic variation, and describe the 
populations from the Gulf of Thailand as a new subspecies. 


Material and methods 

Tissue samples comprising skin and a small quantity of underlying muscle were collected from seven Batagur 
caught in the Sre Ambel River system in Cambodia’s south-west. Samples were taken as a small slither from 
the trailing edge of the outer toe of the rear right foot. The samples were preserved in 75% ethanol for 
shipment and storage at -20°C. These samples are housed in the tissue collection of the Institute for Applied 
Ecology, University of Canberra (voucher numbers AA1001-AA1004, AA1014-AA1016). For one of these 
samples, three mitochondrial and three nuclear genomic fragments were sequenced on an ABI 3130 Genetic 
Analyzer, following the procedures and using the primers described in Le et al. (2007) and Fritz et al. (2008). 
Mitochondrial sequences were for the cytochrome b gene (cyt b) plus the adjacent portion of the tRNA-Thr 
gene (1169 bp), and the partial 12S rRNA (392 bp) and 16S rRNA genes (535 bp). Nuclear sequences 
corresponded to the partial C-mos (590 bp), Ragl (647 bp), and Rag2 genes (672 bp). Sequences of the 12S, 
16S, C-mos, Ragl, and Rag2 fragments were also produced for B. baska, using a sample from the tissue 
collection of the Museum of Zoology Dresden (MTD T 3088, Sunderbans, West Bengal, India), for which the 
cyt b gene was previously sequenced (Praschag et al. 2007). For the other six Cambodian samples and 
additional samples of B. affinis (MTD T 5673-5675) and B. baska (MTD T 5672), only the 
phylogeographically highly informative mtDNA fragment including the cyt b gene was generated. The new 
samples of B. affinis were from confiscated terrapins of Indonesian origin; the new B. baska sample was 
obtained from a market specimen from Mongla, Bangladesh. GenBank accession numbers of the 21 new 
sequences are: FN256231-FN256247 and FN313567-FN313570. 

Sequences from previous studies (Le et al. 2007; Praschag et al. 2007, 2008) were downloaded from 
GenBank and aligned with our data in MEGA 4.0 (Tamura et al. 2007). To elucidate phylogenetic relationships, 
the six-gene data set of Le et al. (2007) was supplemented with sequences of B. baska and the Cambodian 
Batagur and their analyses were repeated (it should be noted that sequences labelled by Le et al. 2007 as ‘R. 
baska' are actually B. affinis). The appropriate substitution model for the concatenated sequence data was 
estimated using MODELTEST 3.04 (Posada & Crandall 1998) and MrMODELTEST (Nylander 2004). According 
to the Akaike information criterion, the best fit model was GTR+I+G with empirical base frequencies: 
7 iA= 0.2943, 7tC=0.2670, 7iG=0.2074, 7iT=0.2313; proportion of invariable sites 1=0.5723; gamma shape 
parameter a=0.5496; rate matrix: R(a)[A-C]=3.5980, R(b)[A-G]=l 1.9253, R(c)[A-T]=4.9027, R(d)[C- 
G]=1.5067, R(e)[C-T]=36.1455, R(f)[G-T]=1.0000. Phylogenetic trees were reconstructed using Maximum 
Parsimony (MP) and Maximum Likelihood (ML) as implemented in PAUP* 4.0b 10 (Swofford 2002) and 
Bayesian inference of phylogeny as implemented in MrBAYES 3.1.2 (Ronquist & Huelsenbeck 2003). Two 
Rhinoclemmys species (Le et al. 2007) were used as outgroups. Both MP and ML analyses were performed in 
a heuristic search with TBR branch swapping option (10^ rearrangements) and best-fit model settings applied 
to the data set in ML. Under MP, gaps were coded as fifth character state; 3155 of 4015 aligned characters 
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(including gaps) were constant in the ingroup sequences; 560 characters were variable and parsimony- 
informative; 300 variable characters were singletons. When the two outgroup sequences were considered, 
3093 sites were constant; 629 were variable and parsimony-informative and 293 variable characters were 
parsimony-uninformative. Clade support for MP trees was estimated using 1000 bootstrap replicates 
(Felsenstein 1985) in a fast heuristic search with all characters unordered and equally weighted and gaps 
treated as fifth character state. ML bootstrap support was obtained by 100 bootstrap replicates in GARLI 0.951 
(Zwickl 2006) with the model parameter raw string defined within a separate GARLI nexus block file. Using 
default settings, branch lengths were optimized with a threshold of 2 x 1 O'* generations per bootstrap replicate. 
Constant -InZ values were reached after a few hundred generations and the allowed minimum value of 
optimization precision was reached after approximately 5500 generations in different replicates. Therefore, 
the parameter <genthreshfortopoterm> was set to 6000 in order to speed up the analysis. Two independent 
search replicates per bootstrap replicate were performed by default settings <searchreps=2> to increase the 
chance of finding the best tree per bootstrap replicate. 

Bayesian analyses (BA) were performed using the Metropolis-coupled Markov chain Monte Carlo 
algorithm with two parallel runs, each with one cold and three heated chains. In a first BA the concatenated 
sequence data set was divided into six partitions corresponding to the three mitochondrial and three nuclear 
fragments and the GTR+I+G model was applied to each partition. The overall rate was allowed to vary 
between partitions by setting the priors <ratepr=variable> and model parameters (gamma shape, proportion of 
invariable sites, etc.) unlinked across partitions, so that for each partition a separate set of parameters was 
estimated. In addition, a second unpartitioned analysis was run. For both analyses, the heating parameter X 
was set to 0.1 to obtain convergence. The chains ran for 10® generations with every 100* generation sampled 
(burn-in=5000). The remaining trees of each analysis were used for generating 50% majority rule consensus 
tree. The posterior probability of any individual clade in such a consensus tree corresponds to the percentage 
of all trees containing that clade, and is by thus a measure for clade frequency and credibility. 

To explore differentiation of mitochondrial haplotypes, parsimony networks were constructed for two data 
sets of cyt b using TCS 1.21 (Clement et al. 2000). The first data set comprised short aDNA sequences (320 
bp) of historical type and museum specimens of B. affmis and B. baska (Praschag et al. 2008) combined with 
the corresponding fragments of sequences generated in the present study and GenBank sequences for B. 
affinis, B. baska, and B. kachuga. The second, 1067-bp-long data set (1038 bp cyt 6 + 29 bp tRNA-Thr) 
included only the GenBank sequences of this minimum length and sequences of the present study. Sequences 
AY434600 {B. affinis) and EU030215 {B. kachuga), which lacked the terminal 29 bp of tRNA-Thr, were 
placed in the last two positions of the alignment. For this second data set, uncorrected p distances were 
calculated in MEGA 4.0. 

External morphology of Cambodian Batagur was compared with museum and live specimens and photos 
of B. baska from the Sundarbans (India, Bangladesh) and Myanmar, with B. affinis from Indonesia 
(confiscated animals) and the west coast of peninsular Thailand and Malaysia, and with B. affinis from the 
east coast of peninsular Malaysia. In addition, one subadult male Irom the questionable locality ‘Indus Delta, 
Sindh’ (Pakistan) from the collection of the Natural History Museum Vienna was examined. For catalogue 
numbers of studied museum specimens, see Praschag et al. (2008), if not given below. Eive B. affinis and 
photos of live terrapins were Irom the Klong Ea-ngu River (Satun Province), Thailand and from the Dungun, 
Perak, and Terengganu Rivers, Malaysia. The Klong Ea-ngu and Perak Rivers are on the west coast, the 
Dungun and Terengganu Rivers on the east coast of the Malay Peninsula. 


Results 

Phylogeny. The topologies of our BA, ME and MP trees were broadly concordant. The BA and ME trees 
were identical, and the single most parsimonious tree (2282 steps; CI=0.521, RI=0.517, RC=0.270) differed 
only in the placement of Siebenrockiella crassicollis. In the MP analysis, this species was basal to all other 
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geoemydids of the ingroup, albeit with weak bootstrap support (<50%), whereas in the BA and ML analyses 
the two Geoemyda species were basal (Fig. 1). Our phylogeny, based on the same three mitochondrial and 
three nuclear genes, had the same general topology as the phylogeny previously reported by Le et al. (2007) 
but with the addition of true B. baska and the Cambodian Batagur. 

With respect to Batagur, our phylogenetic analyses demonstrate, with high statistical support, that the 
riverine B. kachuga is basal to a terminal clade comprising the estuarine B. baska, B. affinis, and the 
Cambodian Batagur. This clade is sister to another clade consisting of the other three Batagur species, B. 
dhongoka + (B. borneoensis + B. trivittata). 


^ estuarine 
I I riverine 
15 both 


0.1 


Rhinoclemmys melanosterna 
Rhinoclemmys nasuta 



FIGURE 1. Bayesian reconstruction of the phylogeny of Batagur and allied geoemydid taxa, based on the expanded 
data set ofLe et al. (2007). Numbers above nodes are posterior probabilities (partitioned analysis). Posterior probabilities 
are identical for unpartitioned analysis, except for the clade comprising B. dhongoka + (B. borneoensis + B. trivittata) 
and the basal clade of all taxa except Geoemyda and Rhinoclemmys (.97 and .81, respectively). Numbers below nodes, 
ML and MP bootstrap values. Vox Pangshura + {Hardella + Batagur) the habitat is coded. 

Haplotype networks and uncorrected p distances. For the 320 bp data set and 90%-95% connection 
limits, the parsimony network analyses yielded three unconnected networks—one for each of Batagur baska, 
B. kachuga, and B. affinis plus the sequences of the Cambodian Batagur (not shown). The haplotype of the 
Cambodian Batagur was separated by seven mutational steps from the more frequent haplotype of B. affinis 
with the latter differing from the rarer haplotype of B. affinis by one step. Within the other two species 
haplotypes differed by a maximum of three steps each. When a connection was enforced, the number of 
mutational steps did not change in each of the resulting three subnets. However, haplotypes of B. affinis were 
joined then via a loop, so that the Cambodian haplotype was alternatively also connected via nine mutational 
steps with the rarer haplotype of B. affinis (Fig. 2a). 

Using the 1067 bp data set, four unconnected networks or single haplotypes were obtained under the 95% 
threshold (one for each of Batagur baska, B. kachuga, B. affinis, and the Cambodian Batagur, not shown). 
Compared to the 320 bp data set, one additional mutational step occurred within the net of B. kachuga, so that 
an ancestral haplotype was separated by two steps from each of its descendants. Under the 90% threshold, the 
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Cambodian haplotype was connected via 18 steps with each of the B. affinis haplotypes, while the haplotype 
of B. baska and the net of B. kachuga remained separate. When a connection was enforced, haplotypes of B. 
baska and B. kachuga differed by a minimum of 50 mutational steps and a maximum of 54 steps; haplotypes 
of B. baska and B. affinis were separated by 44 to 47 steps; and haplotypes of B. affinis and B. kachuga by a 
minimum of 54 and a maximum of 61 steps. The Cambodian haplotype differed from the closest other 
haplotypes, the two haplotypes oiB. affinis by 18 to 19 steps (Fig. 2b). This pattern is also echoed by the 
uncorrected p distances (Table 1). The Cambodian Batagur sequences differed on average from those of B. 
affinis by 1.692% and from B. baska and B. kachuga by 4.592% and 5.635%, respectively. The latter values 
resemble the uncorrected p distances occurring between B. affinis, B. baska and B. kachuga (4.220%- 
5.085%). 



FIGURE 2. Parsimony networks for mitochondrial haplotypes oiBatagur affinis, B. baska, B. kachuga, and the 
Cambodian Batagur (connection enforced). Symbol size corresponds to haplotype frequency; missing node haplotypes 
black. Lines joining haplotypes, one mutational step except otherwise indicated, (a) Network based on a 320-bp-long 
alignment of cyt b. Haplotypes and their frequencies (see Appendix): B. affinis - A1 («=9), A2 (n=l); B. baska - B1 
(«=5), B2 (n=l), B3 («=1); B. kachuga - K1 (n=3), K2 (n=l), K3 («=1); Cambodian Batagur - C (n=7). Haplotypes 
A1 and B1 include the lectotype of Tetraonyx affinis Cantor, 1847 and topotypic specimens of Emys baska Gray, 1830, 
respectively (Praschag et al. 2008). Haplotypes K1 and K2 are from topotypic specimens of Emys kachuga Gray, 1831. 
(b) Network based on a 1067-bp-long alignment of cyt b. Haplotypes and their frequencies: B. affinis - Al («=8), A2 
(«=l);fi. baska-'&l («=4); i?. kachuga - (n=3), K2 (n=l), K3 («=1); Cambodian jSatogwr-C (n=7). 

External morphology and natural history. Morphologically, Batagur baska is highly distinct from B. 
affinis (Praschag et al. 2007, 2008). Batagur affinis from the east and west coasts of the Malay Peninsula are 
known also to differ morphologically (Moll 1980; pers. observ.; see also Figs 3-4; Table 2); our sequence data 
of B. affinis all correspond to the west coast form that seems also to occur in Sumatra (cf. the identical 
mtDNA sequences of confiscated Indonesian terrapins from Praschag et al. 2007 and this study; Sumatra is 
the only part of Indonesia where B. affinis occurs). Cambodian Batagur resemble B. affinis from the east 
coast, but males slightly differ in coloration, so that most adult males can be reliably distinguished. 
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FIGURE 3. (a) Batagur baska, male, Sundarbans, Bangladesh - photo: S.M.A. Rashid; (b) B. baska, semiadult female 
(the pointed, upturned snout develops only with increasing age), Sundarbans, Bangladesh - photo: R Praschag; (c) west 
coast form of B. affinis, male, Klong La-ngu River, Satun Province, Thailand - photo: B. Horne; (d) west coast form of 
B. affinis, female, Perak River, Malaysia - photo: E.O. Moll; (e) east coast form of B. affinis, male, Dungun River, 
Malaysia - photo: E.H. Chan; (f) east coast form oiB. affinis, female, Terengganu River, Malaysia - photo: E.O. Moll; 
(g) Cambodian Batagur male, Sre Ambel River system, Cambodia - photo: R. Holloway; (h) Cambodian Batagur 
female, Sre Ambel River system, Cambodia - photo: B. Horne. Note differences in head shape, soft part and iris 
coloration. 
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FIGURE 4. Hatchlings of Batagur affinis, (a) west coast form, Perak River, Malaysia; (b) east coast form, Terengganu 
River, Malaysia - photos: E.O. Moll. Note yellow marginal scutes and silvery blotches in temporal and parietal region in 
the east coast hatchling. 

TABLE 1. Uncorrected p distances (percentages) within and between Batagur affinis, B. baska, B. kachuga, and the 
Cambodian Batagur based on a 1067-bp-long alignment of cyt h. n, number of haplotypes. Below diagonal, average 
between-group differences; on diagonal, average within-group differences in bold. In brackets, standard error estimates 
(500 bootstrap replicates). 



n 

affinis 

baska 

kachuga 

Cambodia 

affinis 

2 

0.021 (0.020) 




baska 

1 

4.220 (0.635) 

0(0) 



kachuga 

3 

5.085 (0.703) 

4.937 (0.667) 

0.153 (0.071) 


Cambodia 

1 

1.692 (0.379) 

4.592 (0.634) 

5.635 (0.720) 

0(0) 


These morphological differences are paralleled by differences in nesting behaviour. Batagur baska nests 
along the sea shore or on sandy islands in the brackish estuaries, while B. affinis and Cambodian terrapins 
swim far upstream for nesting to reach beaches that are often located well above tidal influence. Some B. 
affinis females may swim over 80 km to these sites (Moll 1980; Ernst et al. 2000; Kalyar et al. 2007). On the 
west coast of peninsular Malaysia, females lay each clutch in a single nest and then make a false nest (‘body 
pit’) nearby. In the Terengganu River population (east coast), by contrast, the females often divide the clutch 
into two or three separate nests (Moll 1980 and pers. comm.). 


Taxonomic conclusions 

Our genetic and morphological data demonstrate the distinctness of the disjunct Cambodian population of 
Batagur. However, the degree of genetic differentiation is clearly lower than that between B. affinis, B. baska, 
and B. kachuga (Fig. 2; Table 1). The Cambodian Batagur is undoubtedly very closely related to B. affinis, 
which we interpret as subspecific rather than interspecific level variation. 

Although slight coloration differences exist between the two (Table 2), the Cambodian Batagur resembles 
morphologically populations of B. affinis from eastern Malaysia (Dungun and Terengganu Rivers) that were 
unfortunately not available for genetic study. We suggest that these populations and the Cambodian Batagur 
represent the same taxon and became isolated owing to human activity. We contend that this taxon was 
distributed originally in estuaries surrounding the Gulf of Thailand. 

The type locality of Penang for B. affinis is along the west coast of peninsular Malaysia which, together 
with mtDNA sequence data of the lectotype (Praschag et al. 2008), provides evidence that the name Tetraonyx 
affinis Cantor, 1847 refers to the west coast taxon and that no name is available for the east coast form, also 
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distributed in Cambodia (cf. Fritz & Flavas 2007; Praschag et al. 2008). Consequently, we deseribe the eastern 
form of B. affmis as a subspecies new to science. 


TABLE 2. External morphology and current distribution ranges of estuarine Batagur based on personal observations and 
data from Anderson (1879), Rashid & Swingland (1997), and Moll (1980 and pers. comm.). 


Batagur baska 


Batagur affinis west coast Batagur affinis east coast Cambodian Batagur 


Male 


Female 


Hatchling 


Range 


Head elongated with 
pointed, upturned snout; 
head black, area around 
nostrils pale bluish, rest of 
head and distal part of neck 
deep black, passing into rich 
crimson on base of neck; iris 
matte greenish yellow; 
forelimbs brilliant rosy 
carmine, hind limbs, tail and 
thighs dull reddish purple. 
Carapace during mating 
season rich brown to 
reddish, in some individuals 
slightly marbled with darker 
lines; plastron having a rosy 
yellow tint 

Head elongated with 
pointed, upturned snout; 
light grey with waxy-blue 
nostrils and yellowish jaws, 
temporal region paler, 
almost white; other soft 
parts grey, never brownish 


Unknown 


Northeastern India from 
Orissa to West Bengal, 
Sundarbans of Bangladesh 
and coastline of Myanmar, 
perhaps northernmost west 
coast of peninsular Thailand 


Head distinctly shorter 
than in B. baska, with 
blunt snout and shorter 
distance from nostril to 
eye; head jet black or very 
dark grey, never 
brownish; iris during 
mating season immaculate 
white. 

Carapace black or dark 
grey, during mating 
season black 


Head short with blunt 
snout; distance from 
nostril to eye very short 
(resembling the riverine 
B. kachuga); head dark 
grey or brownish, jaws 
dirty yellow or light 
brown; other soft parts 
dark grey or brownish 

Overall dark grey, without 
light temporal blotches 


Southwest coast of 
peninsular Thailand, west 
coast of peninsular 
Malaysia, Sumatra 
(Indonesia) 


Head elongated with 
pointed, upturned snout and 
in comparison to west coast 
males longer distance 
between nostril and eye; 
head and soft parts chocolate 
brown to almost black, light 
coloured individuals with 
brown and never grey skin; 
edges of mouth orange; iris 
golden or bright yellow. 
Carapace dark brown to 
black 


Head elongated with 
pointed, upturned snout like 
in B. baska', head greyish to 
brownish with whitish grey 
to silvery blotches in 
temporal and parietal region 
and brown jaws 


Overall brown, with whitish 
grey to silvery temporal and 
parietal blotches and distally 
yellow marginal scutes 

East coast of peninsular 
Malaysia and of peninsular 
Thailand 


Closely resembling 5. affinis 
from the east coast of 
Malaysia. However, head 
often rusty brown to reddish, 
distal portion of neck turns 
proximally into greyish; 
limbs greyish as well. 
Carapace dark grey 


Look exactly like females 
from east coast of Malaysia 


Look exactly like hatchlings 
from east coast of Malaysia 


Sre Ambel River system, 
Cambodia 


Batagur affinis edwardmolli subsp. nov. 

Holotype. Natural History Museum Vienna, NMW 38903, juvenile in aleohol (hatehed and died in eaptivity), 
Sre Ambel River system, Koh Kong Provinee, Cambodia; don. Head Start Centre Sre Ambel, July 2009. 

Paratypes. Museum of Zoology Dresden, MTD 47538, juvenile in aleohol, same data as holotype. Field 
Museum of Natural History, Chieago, FMNH 224093 (ex EOM 2398; speeimen figured on Plate IIIA in Moll 
1980), broken shell of adult female, Terengganu River, Malaysia; leg. Edward O. Moll, 7 July 1976. 

Etymology. The new subspeeies is named in reeognition of Professor Edward O. Moll, one of the 
foremost experts on river turtles, who substantially eontributed to the knowledge of Batagur affmis and its 
natural history. 

Diagnosis. Adults differ from nominotypieal subspeeies of Batagur affinis by their distinetly more 
elongated head with upturned snout; males with ehoeolate brown to almost blaek head (east eoast of 
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peninsular Malaysia) or sometimes rusty brown to reddish head (Sre Ambel River system, Cambodia), edges 
of mouth orange; iris golden or bright yellow. Females and juveniles with eonspicuous whitish grey to silvery 
blotehes in temporal and parietal region; hatehlings with distally yellow marginal seutes. For eorresponding 
eharaeters of B. a. qffinis, see Table 2. 

Description of holotype. Speeimen slightly maeerated; some epidermal seutes detaehed from shell. 
Carapace roundish when viewed from above, with weakly serrated central and posterior marginal scutes; 
medial keel distinct, with posteriorly directed, slightly pointed spines. Plastron anteriorly truncated, 
posteriorly with anal notch. Straight line carapace length approximately 86 mm, carapace width 84 mm; 
medial plastron length 74 mm, maximum plastron length (to tips of anal scutes) 78 mm. 

Range: East coast of peninsular Malaysia and adjacent Thailand; Sre Ambel River system, Cambodia 
(Fig. 5). 



FIGURE 5. Historical distribution of Batagur affinis affinis, B. a. edwardmolli, and B. baska (modified from Praschag et 
al. 2008). Note that the species are extirpated in most of their former ranges. 
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Discussion 


Our phylogenetic analyses suggest repeated switches of habitat during the evolution of the extant Batagur 
species. The vast majority of the some 60 geoemydid species (Fritz & Flavas 2007) are freshwater terrapins, 
while B. affinis, B. baska, and B. borneoensis are the only species living in tidal, brackish areas of the 
estuaries of medium and large rivers (Ernst et al. 2000). In addition, the nearly extinct B. trivittata once used 
this habitat, but also lives and nests until today far upstream in the inland of Myanmar (Maxwell 1911; 
Kuchling et al. 2006). When the natural history of the closely related genera Hardella and Pangshura is 
considered, the most parsimonious view has the ancestor of Pangshura, Plardella, and Batagur as a riverine 
freshwater terrapin (Fig. 1). If this is true, the ancestor of the Qsdani Batagur taxa may well have been riverine, 
in which case an estuarine mode of life has been independently acquired in each subclade of Batagur by B. 
affinis, B. baska, B. borneoensis, and B. trivittata. Consequently, the riverine mode of life of the other species 
would represent the ancestral ecological adaptation. Alternatively, the ancestral Batagur could have been an 
estuarine terrapin and B. kachuga, B. dhongoka, and B. trivittata may have returned independently to the 
riverine freshwater habitat. 

yiosX Batagur species exhibit conspicuous sexual and seasonal dichromatism (Anderson 1879; Moll 1980; 
Moll et al. 1981; Praschag et al. 2007)—an extremely rare trait among chelonians. According to Klingelhoffer 
& Mertens (1944) and Moll et al. (1981), the bright coloration of males is associated with reproduction. A 
function as premating isolation mechanism seems plausible given the capability of many chelonians to 
hybridize (cf Fritz et al. 2008) and that, where the similar-sized Batagur species occur syntopically, they 
differ considerably in breeding coloration. For instance, B. affinis often occurs together with B. borneoensis. 
During the mating season, the coloration of males of the two species is highly distinct. Breeding males of B. 
borneoensis have a light coloured carapace with three broad black longitudinal stripes; the skin of the head 
gets bright white with a conspicuous black-edged scarlet stripe running between the dark eyes from the nose 
to the occiput (Moll et al. 1981). On the west coast of the Malay Peninsula, males of B. a. affinis have instead 
a black carapace and a jet black head with immaculate white eyes; on the east coast, the head of males of 5. a. 
edwardmolli gets chocolate brown to black, with orange edges of the mouth and a golden iris. These 
contrasting differences are suggestive of character displacement (Moll et al. 1981). The coloration differences 
of the two B. affinis subspecies and of B. baska (Table 2) could also act as isolating mechanism between these 
taxa. While the two B. affinis subspecies are entirely allopatrically distributed, there is at least the possibility 
of a historical sympatric occurrence of B. baska and B. a. affinis near the border region of Myanmar and 
peninsular Thailand. 

Also with respect to nesting behaviour, character displacement seems to occur between different Batagur 
species. Batagur affinis, which occurs sympatrically with the estuarine B. borneoensis, nests far upstream, 
whereas B. borneoensis nests on the same sea beaches as marine turtles (Dunson & Moll 1980). By contrast, 
B. baska, occurring farther northwest where B. borneoensis is absent, uses sea beaches for nesting (Ernst et al. 
2000). However, in the Ayeyarwady Delta, large numbers of B. baska and B. trivittata nested together until 
about 100 years ago (Maxwell 1911; Kuchling et al. 2006). 
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Appendix. Mitochondrial cyt h haplotypes oiBatagur taxa from Fig. 2 (GenBank accession numbers; short sequences of 
320 bp asterisked). For new sequences the collection vouchers are given in brackets (AA - Institute for Applied Ecology, 
University of Canberra; MTD T - Museum of Zoology Dresden, Tissue Collection). 

Batagur affinis—Al («=9): AM691750-AM691754, FN313568-FN313570 (MTD T 5673-5675), AM922509*; A2 
{n=\y. AY434600. 

Batagur baska— («=5): AM495267-AM495269, FN313567 (MTD T 5672), AM922507*; B2 (« = 1): AM922510*; 
B3 (« = 1): AM922508*. 

Batagur kachuga—YA («=3): AM495284-AM495285, AM495287; K2 (n=l): AM495286; K3 («=1): EU030215. 
Cambodian fiatagur—C (n=7): FN256231-FN256237 (AA1001-AA1004, AA1014-AA1016). 
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